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ABSTRACT 


A three-dimensional, primitive equation model by Anthes 
me Warner (1978) was used to predict five typhoon cases. 

The five-layer model on a 40X40 staggered grid with 120 km 
resolution included a moisture cycle, sensible and latent 

heat flux at the earth's surface, and a bulk parameteriza- 

tion of the planetary boundary layer. The model is initialized 
uSing operational wind fields and two forms of a wind bogus. 

Due to the lack of a representative moisture analysis 
mieene Vicinity Of tropical cyclones, simulated moisture 
fields were used to initialize the model. Initial experiments 
conducted with these fields produced widespread convection 
and heating which developed circulations in areas well removed 
from the actual storm. The associated modifications to the 
Steering flow, and the overly intense storm circulations 
resulted in premature recurvature. Use of a second moisture 
bogus, which provided less available moisture (especially at 
upper levels), reduced the amount of convection over the 
entire grid and the effect on the steering current. 

A second wind bogus based on a scheme described by Hover- 
male (1976) was also tested. This technique resulted in a 
Smebor end less intenSe initial storm, which also resulted 
in considerably less convection. The tracks forecast by 
the different wind and moisture fields were compared with 
forecasts of the 60 km resolution Madala and Hodur model 


uSing the same initial data. 
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i) LNPRODUCTION 


The major objective of this set of experiments was 
to test a general mesoscale model using operationally- 
analyzed data to determine effectiveness in forecasting 
tropical cyclone development and movement compared to some 
emer tropical cyclone models currently being tested. The 
Penn State model was selected because it contained boundary 
Byer Physics, a Convective moisture cycle, and other 
Options which could be easily included in (or excluded 
from) any forecast. It was then possible to determine the 
effects and relative importance of the various parameteriza- 
peems by Selectively incorporating different combinations. 
This model was also chosen for its general code which allowed 
changes to be made easily to the resolution and size of 
the vertical and horizontal domains. These changes were 
necessitated by time and computer availability constraints. 

Bering the initial forecasts, 1t was discovered that 
the initial wind and moisture conditions played a very 
important role in determining the development and movement 
of tropical cyclones. An investigation was then conducted 
to determine the extent of the effects of the initial wind 
and moisture conditions on the model forecasts. This paper 
presents the results of these experiments and compares fore- 
casts made using the Penn State model with forecasts made 


using two other tropical cyclone models. 
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Poe sen le LLONVOr “tHe MODEL 


The experiments used a general, three-dimensional, 
primitive equation model designed to forecast flows with 
characteristic wavelengths of 10-2500 km. Provisions 
of the model included a variable terrain, a moisture cycle, 
latent heating, sensible heat flux at the earth's surface 
and boundary layer physics. The moisture cycle modeled 
both stable (nonconvective) precipitation with grid-scale 
relative humidities in excess of 100%, and unstable (convec- 
eye) precipitation with grid-scale relative humidities 
less than 100%. Details, such as the finite difference 
equations and staggered grid lattice, were retained intact 
as described by Anthes and Warner (19 78). 

iiemeareicular version of the model used in thas set 
of experiments was similar to one described by Anthes 
(1978). Model inputs consisted of pressure weighted wind 
components (p*u,Psv), relative humidity and temperature at 
five levels, surface pressure, sea surface temperature and 
terrain. Predicted parameters were wind components (u,v), 
Surface pressure (D.), aa eneae ie elaleiub oben ame Gulee is ebbas) —  Sraiel 
potential temperature (9). Anthes's(1978) version was 
reduced from 7 levels to 6 levels at which ¢ was defined 
fo oer Ons, 0.5,0-7,0.9 and 1.0) with all other variables 
(V,T,O,W, and q) defined halfway between these levels at 


oe ormls)054,0.6,0.8 and 0.95. Sigma level 0.0 corresponded 


e3 








to a pressure of 0 millibars. To reduce computer time, 





a 40x40 grid in Mercator coordinates with 120 km resolution 
was selected rather than the 50x50 grid with 60 km resolution 
used by Anthes (1978). Other simplifications were made to 
the model for economic reasons as described by Anthes (1978). 
A spatially and temporally constant vertical distribution 
| of latent heating associated with cumulus convection was 
Specified in lieu of a one-dimensional cloud model (Anthes, 
mee}, and a Simplified version of Deardorff's (1972) bulk 
parameterization of the planetary boundary layer (PBL) was 
substituted for the high resolution PBL (Anthes and Warner, 
1978). This scheme assumed that the surface wind, Vel, 
was equal to the wind at the lowest level (0 = 0.95). 
All experiments described here used an ocean terrain with 
Pembina, and aS Such, a constant value of eS LOLS was 
used for the drag coefficient, Ch: 
All layers were checked for supersaturation and ¢ and 
q were adjusted to maintain relative humidities less than 
100% by condensing a fraction of the moisture excess and 
adding the associated latent heat to the temperature ten- 
dency. The original model assumed a constant fraction of 
-31 for all levels. This factor was modified using an 
approximation suggested by Haltiner (personal communication). 
A new constant was calculated for each level usSing the 
equation 
Shee 


Aq ae ih noe ae (1) 
[l+L se G3 Ro ae | 
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where typical values for q and T were obtained from 
Jordan's (1958) mean September tropical sounding and 
Sheets' (1969) mean hurricane sounding. The modified 
values which were equal to the reciprocal of the denomina- 
mor im equation (1) are listed in Table l. 

Initial conditions were derived from the Fleet Numerical 
Weather Central (FNWC) Global Band Analysis. Initiali- 
zation was accomplished using a scheme developed by Hodur 
(personal communication) for uSe with the Madala and Hodur 
77) split Semi-implicit tropical cyclone model. This 
scheme used a reverse balancing technique to derive interior 
geopotential values on Sigma levels from the nondivergent 
wind components using a nonlinear balance equation. Balancing 
was conducted over a 51x51 grid with 120 km resolution, 


from which the 40x40 grid was interpolated for input to 


the model. 
Table 1. Fractions used for condensing excess 
moisture to obtain 100% relative 
Nuntatevewe values fer T and gq were 
taken from representative tropical 
soundings. Values of q_ were calcu- 
lated using the Clausius-Clapeyron 
equation. 
Level (3) | (°K) qe (gkg +) Fraction 
O15 ZO. 71 6 0.0433 Or o.86 
0.40 IBS 8) VAS ey 0) oS) 2) Oreo Lo 
02.60 25 t0 6 1. OSoe O74 22 
One 0 Zao. 36 12.0868 O3:35 
Oo > 295.76 Hs=s00e8 Cr 30 


digs. 





Initial wind conditions were composed of both analyzed 
and idealized values. Following Hodur (personal . 
tion), a symmetrical wind bogus (referred to as V-bogus 1) 
was obtained within a 360 km radius of the storm center by 
adding idealized values to the vorticity field (which was 
derived from the analyzed winds) prior to calculation of 
the stream function and nondivergent winds. Bogus vorticity 


values were calculated for the surface to 700 mb from 





where are was the magnitude of the wind at the radius 

of maximum wind (TO). The vorticity maximum at 100 mb was 
set equal to Oss G02 sl, Withmt@e: VOreEICiEylat inter 
mediate levels varied linearly between 700 mb and 100 mb. 


Horizontally, the vorticity was decreased in accordance 


wich 


where r was the radial distance (km) from the storm center. 
Values for ae were taken from Typhoon Reports and typhoon 
warning messages and the radius of maximum wind (RMW), to: 
was set equal to 60 km. 


Due to the lack of representative moisture analyses in 


Siemvictartcey Of cropical cyclones, idealized moisture fields 
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(referred to as Rh-bogus 1) were used to initialize the 
model. These fields were produced by specifying a uniform 
field of large scale relative humidities which decreased 
with height. Further bogusing of the storm was accomplished 
by centering a five grid point square of higher relative 
humidities, which also decreased with height, over the 
Storm. In this manner, the storm started with more mois- 
ture than the surrounding environment. Large scale and 


storm relative humidities are given in Table 2. 


Table 2. Relative humidities used in Rh-bogus 1 and Rh-bogus 2 


Level (c) SEOmm Large Scale 
Oe5 2 50 250 
0.40 . 3 oO 
0.60 90 aah 
Ome 0 0) aod 
O95 S/S wo 


Relative humidities used in Rh-bogus 2. 


Level (3c) Seon Large Scale 
low igilteysi 

ORR es 250 -40 Ae 

0.40 Rho ke Eo 0 -46 

Or. 6:0 Evie Aaeiey| py al 

O80 , aoe 16) & 

Opa95 - 80 4 S18) 50 
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Pie EXPERIMENTS 


The basic experiment consisted of 48-hour forecasts 
for the following five cases using all of the model's 
Options (with the exception of terrain); 

i Typhoon OLIVE 78042112 (year,month,day,12 GMT) 

ole Typhoon OLIVE 78042400 

BEt. Typhoon KIM a EO oT 

LEAVE Typhoon PAMELA 76051900 

Wie Typhoon PAMELA 76052200. 
Tracks of these storms were then compared to the tracks 
obtained from forecasts made with the Madala and Hodur 
(1977) split semi-implicit (SSI) model, the FNWC Opera- 
tional Tropical Cyclone (Channel) model (Mihok and Hinsman, 
i7/), ana the best track analyses. Appendix A contains a 
comparison of these two models with the Penn State model. 
Further experiments were then conducted to test the model's 
SeGioteryilty tO the iniltlal wind and moisture conditions and 


to the convective parameterization. 


A. Seowmms OF THE FIRST COMMON EXPERIMENT 

All five cases experienced strong recurvature. See 
Preseme7 through 21 for positions relative to the best track. 
Examination of the wind and moisture fields revealed several 
characteristics which were common to all cases. Every storm 
developed strong frictional convergence with large cross- 


lsobaric flows at low levels, even in areas of high wind 


18 





velocities. Although the maximum storm intensities varied, 
all five cases developed a southerly flow at middle levels 
just to the east of the storm. An example of this circu- 
lation, which was especially prevalent at 600 mb for Typhoon 
KIM, can be seen in Fig. 1. A tendency to produce a large 
fm@eteyclonic circulation at 150 mb with strong westerlies 
located along the northern branch of the gyre also existed 
moet cases. An example of this circulation is given in 
Bee 2. This circulation, combined with the strong middle 
Yevel southerly flow, probably contributed significantly 
BOnbecurving the storms. 

Another dominant characteristic common to all cases, 
which contributed to the recurvature effect, was the produc- 
Seton Or widespread convection, with particularly strong 
Semvective activity Occurring after 12 to 24 hours in the 
Merpedeast Corner and along the northern boundary. All 
cases experienced excessive vertical velocities accompanied 
by strong latent heating and precipitation (both convective 
and nonconvective) in these areas, which resulted in the 
formation of localized areas of lower pressures. Continued 
development of a trough in the northeast corner distorted 
the entire surface pressure pattern to such an extent that 
the storm was steered toward the corner, producing even 
greater development in this region and more rapid recurva- 
ture of the storm. An example of this sequence can be 
seen in Figs. 3a through 3d. Forecasts for both cases of 


Typhoon OLIVE also developed severe convection in the 


iO 








southwest corner and an intense circulation was produced 
Movoemec4 and 36 hours. The first OLIVE case also 
Geveloped an additional circulation which had existed 
initially to the west of the storm. An example of the 
areas which underwent development due to convective activity 
Gam bse seen in Figs. 4 through 6. The development of these 
circulations near the boundaries prompted an investigation 
into the moisture fields, latent heat parameterization 
scheme and the boundary conditions. These topics are 
examined in the following sections. The general procedure 
was to select a typical case and examine successively the 
SeeeeceoOr each modification on that storm. After adopting 
the modifications based on the one case, the entire set of 


mye Cases was re-run. 


BeeeeMOdSTURE EXPERIMENTS 

A new moisture initialization was introduced in an attempt 
to reduce the amount of convection and high relative humidi- 
ties at upper levels in regions not associated with the 
storm. The second technique (referred to as Rh-bogus 2) 
also consisted of an area of higher relative humidities 
Superimposed on a uniform field of lower values. However, 
the second moisture bogus covered an area with a radius of 
3 grid points (360 km), which was the same area covered 
by the wind bogus. The large scale values, which were 
based on Jordan's (1958) mean September tropical sounding, 


varied only with height, while the storm relative humidities, 


20 





taken from Sheet's (1969) mean hurricane sounding, depended 
on the storm surface pressure as well as height. Values 
used in this scheme are listed in Table 2. The high range 
of relative humidities was used for storm surface pressures 
less than 995 mb and the low set for pressures greater than 
1014 mb, with a linear interpolation between those sets 

Hoe intermediate storm surface pressures. Note that the 
major difference in Rh-bogus 2 was the reduction of all 
relative humidities, especially at upper levels. Additional 
48-hour forecasts were then conducted for both cases of 
OLIVE using the new moisture bogus. OLIVE was chosen 
because it had the most widespread convection of all the 


Casco. 


Seeeseoturs OF RH=BOGUS 2 

The results of the low moisture bogus are typified 
in Figs. 7 through 9. Convective activity was reduced over 
the entire grid with the greatest effect noticeable in the 
reduced development of those areas not associated with the 
Storm. The total convective rain produced over the grid 
mas redueed from 382 cm to 198 cm, while the total non- 
convective precipitation was reduced from 26 cm to 7 cm 
in the first 12 hours. More significant, however, was the 
reduction of the maximum convective precipitation associated 
with the second cyclone (located to the west of OLIVE) 
from 12.5 cm to 1.4 cm. The maximum convective precipitation 


associated with OLIVE was only reduced from 5.9 cm to 4.5 cm. 


val 





The associated convective heating rate at 12 hours was 
Heweread trom 280 °C/day to 60 °C/day over the second 
cyclone, while the heating over OLIVE was reduced from 
160 °C/day to 110 °C/day. Similar results were obtained 
im Ene Other areas, particularly in the northeast corner, 
although the development of the trough in this region was 
not prevented. 

@ae net effects of the low moisture bogus were filling 
the storm by about 2 mb and decreasing the winds by less 
than 10%, while significantly reducing the deepening over 
the large scale areas. The general characteristics of the 
flow were not altered, nor were the forecast tracks signi- 


Picancly altered. 


D- NO HEAT EXPERIMENT 

A 24-hour forecast was then made using the high humidity 
bogus, but without including latent heating effects in the 
temperature prediction equation. Typhoon KIM was chosen for 
this experiment because it had the best predicted track of 
the five cases using the SSI model and had been used with 
MieewitnOwt the wind bogus (Hodur, personal communication). 
Deletion of latent heating effects (due both to convective 
and nonconvective processes) resulted in rapid filling of 


the storm and only marginal improvement of the storm track. 


HeeereiGie LO THE MODEL 
Before continuing with further experiments, several 


detailed changes were made to some of the model-defined 
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constants to make the 5-layer version more closely resemble 
the original 6-layer model described by Anthes (1978). 

eee ereical distribution Of convective heating, Nic), 
defined by Anthes (1977), was altered in the lowest two 
layers such that the same amount of heating was provided, 
but over different layer thicknesses. Because this function 
was normalized, the changes produced in these two layers 
WeeeuwcGtouributed throughout the vertical. The corrected 
meoelessOr N(o), which provided the same distribution of 
convective heating described by Anthes (1978) in the 6- 
layer model, are compared in Table 3. A 24-hour forecast 
using the high moisture bogus and the revised heating 
distribution produced a 10% reduction in convection and 

its associated latent heating and precipitation. The storm 
center filled 1 mb. 

Two other changes were made to the convective parameteri- 
Baeron. The vertical eddy fluxes, 6 (q.-d) /80, which 
determined oq/3t in the cloud, were altered to make the 
vertical distribution of cloud moisture more representa- 
tive of the distribution outlined by Anthes (1977). A 
Pe ouestorecast produced less convection, precipitation, 
and heating, and increased the surface pressure OE Che Storm 
by 1 mb. Convection over the entire grid was further 
reduced in a 12-hour forecast by increasing the critical 
required 


g HO ms > 


value of the integrated water vapor convergence, M 
"i 


t! 
6 


for convection to occur from 3.0x10 ° to 3.0x10- 


where 
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Level 


0 


300 


400 


460 


580 


600 


Cow 


800 


820 


940 


250 


Table 3. 


(mb ) 


Revised values of the vertical 
distribution of latent heating, 
Ngee) (New values are given 

for experiment designed to shift 


heating maximum to higher levels.) 


6-layer 


24 


5-layer 
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New 
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i 26 
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The number of grid points experiencing convection for any 
One time step was reduced significantly, with the maximum 
number of grid points (at 12-hours) lowered from 797 to 
Only 222. The total convective precipitation over the grid 
was reduced from 414 cm to 263 cm. The changes in convec- 
tive precipitation and heating produced by the revisions 
discussed here, can be seen in Figs. 10a through 10d. 
Although the motivation for the above changes was to reduce 
what appeared to be excessive latent heating in the model, 
Weme tliat Ene Circulation in the northeast corner still 
experiences strong convergence and excessive precipitation. 
Two experiments designed to test the model's sensitivity 
to the convective parameterization were conducted using Rh- 
Begadse2 £Or initial conditions. A 24-hour forecast was 
made using a vertical distribution of latent heating which 
Shifted the maximum heating to higher levels. This was 
done so that the model would develop more cyclonic circula- 
tion at upper levels and possibly produce a greater westward 
movement. The new values of N(c) are listed in Table 3. 
Additional upper level heating produced less precipitation 
in the northeast corner but failed to alter the forecast 
Significantly in 24-hours. Another 24-hour forecast was 


made with the integrated water vapor convergence, Mey 
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Ga¥culated for only the two lewest layers of the model. 

The objective of this experiment was to remove the effects 
of the large upper level moisture convergence and, thereby, 
reduce the precipitation in the northeast corner while not 
Seriously affecting the storm. Confinement of the moisture 
convergence calculation to low levels reduced convection and 
development of a trough in the northeast corner while 
intensifying convection associated with the storm after 

1i2 hours. However, after 24 hours, the trough in the 
northeast corner developed rapidly and produced results 
Similar to those already encountered. These last two 

model revisions were not incorporated into the model for 
future experiments. 

Frictional effects were also modified for changes made 
in 60 and p* which had resulted from changing the pressure 
at the top of the model from 200 mb to O mb, and decreasing 
the thickness of the surface layer from 150 mb to 100 mb. 
These modifications resulted in a 4% decrease in the effec- 
tive drag coefficient. A 12-hour forecast made using the 
reduced friction deepened the storm surface pressure by 
1 mb and increased the amount of convection occurring over 
the storm, while decreasing the total number of grid points 


experiencing convection. 


Pee SORE ZONTAL DIFFUSION EXPERIMENT 
An attempt was made to reduce the strong horizontal 


shears which existed along the boundaries. An example of 
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eae chear can be seen in Fig. 2 where the interior upper 
evel circulation differs drastically from the conditions 
along the western and southern boundaries. Since this 
Sm@earaceeristic also Geecurred when the model was run with 
extratropical data (see Anthes, 1978, Figs. 23c through 
23e), it was decided to increase the horizontal diffusion. 


Mhe horizontal eddy viscosity, K was multiplied by an 


H’ 
Piplerrercation factor which increased from a value of 1.0 
near the center to a maximum near the boundaries, where 
42 -l1 


Ris Mmanmerbestrictea to be less than 30x10 ms 2 Leah e thera 


runs with the model were made using twice the calculated 

Ko and a maximum value of 60x 10%m-s -. The increased 
diffusion decreased the shears along the boundaries by less 
than 10% and therefore did not significantly reduce the 
convergence which resulted in convection along the northern 
and eastern boundaries. 

Another factor which may have contributed to over- 
development of the northeast corner was the lack of moisture 
diffusion. The winds, especially at upper levels, tended 
mamacdvece MOLSture into that corner. Without diffusion, 
lt may have been possible that moisture simply accumulated 
Deen cerner Producing more convection. This effect can 
Sewseem Dy Comparing Fig. 3c with Fig. 11 which is the 24- 
hour forecast without latent heating. In the latter, although 


convection occurred, the heating erfects were not added to 


the tendency and the circulation near the corner did not 
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develop. As a result, the surface pressure pattern was 


mot distorted. 


G. WIND BOGUS EXPERIMENTS 

As demonstrated in Fig. 3a, the initial wind bogus 
(V-bogus 1) produced a storm which was very large. Although 
bogused values were entered in the vorticity field within 
Suemeadius Of 3 grid points, the cyclonic circulation and 
the associated pressure fields derived from the balance 
equation extended out to a radius of 8-10 grid points. 
It was also desirable to include an inner region of cyclonic 
outflow at upper levels. A second wind bogus (V-bogus 2) 
was calculated uSing the same conditions from the surface 


to 700 mb as V-bogus 1, while at 400 mb, the vorticity was 


given by 
= = rc ah ae 
oy Sax Te ome SS bee sa (4) 
Et = (xr = 180) 
oe = —£ + 130 i Teemu ts os. 60 KM, (6) 
= 0 Ge sO OK, (7) 
where the vorticity maximum, Cmax Was computed in the same 


manner as described earlier for V-bogus 1. At 100 mb, the 


vorticity was calculated using one-half the calculated 
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values for 400 mb. The vorticity at other levels (500, 
300,250,200, and 150 mb) was varied linearly between 700 
and 400 mb and between 400 and 100 mb. This bogus produced 
essentially the same results as V-bogus 1, probably because 
the circulation in the lowest levels of the storm was still 
very large. 

Therefore, a third wind bogus (V-bogus 3) was calculated 
in an attempt to reproduce a wind bogus described by 
Hovermale et al (1976). The tangential wind, Va, was 


calculated from 


ale me oe 
Vg(m s ) = E97 Ou (ee 6.00sem, (8) 
and 
ail _ (r-60 km) 
Vy (m Se) es a Ti) es 10 
60 < ry (9) 


for the lowest level, while the upper level wind was 


derived using 


Soll pa es Lo 
V4 (m s ) = 20 Oa ee oUmikn 7 (10) 
v,(ms *) = 12{1 a= oC ease eG0) km, 


(arin 


2g 





and 


1) (xr-60 km) 


Se -) 660 < 2, (12) 


were r extended out to a radius of 6.4 grid points (768 





Mme The vorticity was then calculated from 
OV V OV 
= Q ee im 
5 = tr * F r 30 ES 


where the last term was taken as zero for a symmetrical 
See@ie Frior tO substituting the bogused vorticity values, 
a Simple Laplacian smoother was passed 30 times over the 
analyzed wind components in the same area covered by the 
bogus. Experiments conducted on both V-bogus 1 and V- 
bogus 3 with and without the smoother did not produce 
pmombetcantiy different results for either bogus. 
The following experiments were then conducted uSing 
Typhoon KIM and Rh-bogus 2; 
ime -SOogus if: | 24=hours 
Mem V-DOgUS 3: 43=hours 
mene V-bogus | with one-half the vorticity 
maximum at the surface as in case I: 
A3=Nnours . 
Differences produced in the initial wind fields by the two 
bogusing techniques can be seen in Fig. 12a through l2d. 
Note that V-bogus 1 produced a very large and intense storm 


compared to the one produced by V-bogus 3. In particular, 
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V-bogus 1 produced large changes in the winds and pressures 
at grid points a considerable distance from the storm. 
Because the 40x40 grid was selected from the interior of 

a 51x51 grid, wind values were changed as much as 50% in 
magnitude and 180° in direction along the boundaries by 
V-bogus 1. Subsequent analysis fields, that were required 
for calculation of the time-varying boundary values, were 
not bogused. The differences in these two fields introduced 
accelerations greater than expected in a 12-hour period 
along the southern and eastern boundaries. These accelera- 


tions were not introduced by V-bogus 3. 


H. RESULTS OF WIND BOGUSING EXPERIMENTS 

Peco l2 hours, V=bpogus 1 produced convection at 100-200 
grid points while V-bogus 3 produced convection at only 20-30 
Sere scolints. Total convective precipitation over the grid 
was lowered from 150 cm to 96 cm and the total stable rain- 
fall was reduced by about 40%. V-bogus 3 dramatically reduced 
the amount of precipitation and heating in the northeast 
corner while concentrating more rain and heating over the 
Storm. The effects of the smaller vortex can be seen in 
Figs. 13 through 15. Note the lack of development of a 
large trough in the northeast corner. The southerly flow 
at middle levels was also reduced. Reduced development of 
these two features resulted in the improved forecast track, 


sieowl ih Fig. 16. 
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The 24-nour forecast developed more convection in the 
northeast corner, however the storm seemed to be suffi- 
Gtentliy far from the corner to not be seriously affected 
By the weak troughing. A large anticyclonic circulation 
developed at 150 mb with strong southerly winds in the 
southwest quadrant. Apparently, these winds contributed 
to the northward movement of the storm. 

Reducing the vorticity of V-bogus 1 to half the 
Original value (Experiment III) resulted in a smaller 
vortex which produced values of convective activity, precipi- 
tation and heating which were between those produced by 
V=bogus 3 and the high intensity V—bogus 1. The track fore- 
case uSing this bogus was also between the tracks produced 
by the other two boguses. This experiment demonstrated 
the effect that variations in the Coriolis parameter can 
have upon excessively large or intense vortices; that 1s, 
larger storms produce larger northward accelerations 


Miovermale, et al., dies ag ae 
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IV. COMPARISON WITH OTHER MODELS 


The effects of V-bogus 1 were compared to the FNWC 
analyzed fields (no storm bogus) for KIM, OLIVE (78042400) 
and PAMELA (76052200) using the SSI model by Hodur (personal 
communication). A 48-hour forecast was also made for 
PAMELA (76052200) using the high humidity moisture bogus 
Peieeegus 1) and V-bogus 1. A summary of the five cases 
available for comparison with experiments uSing the Penn 


meote Model 1S given in Table 4. 


Table 4. Summary of experiments using the 
various bogusing techniques. 


Seomm Sol Penn State 
OLIVE (78042112) Analyzed winds, dry V-bogus 1/Rh-bogus 1 


V-bogus 3/Rh-bogus 2 


OLIVE (78042400) Analyzed winds, dry Same as above 
V—-DOgus sly sanry 


KIM Analyzed winds, dry Same as above 
Ve DOgusml, = cm 


PAMELA (76051900) Analyzed winds, dry same as above 
PAMELA (76052200) Analyzed winds, dry Same as above 


V = DOG Si meray: 
Me-bogus 1, wet 


Forecast tracks obtained from the common experiments are 


Seomm in Figs. 17 through 21. V=bogus 1 produced more 


recurvature in all cases tested with the SSI model, although 
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not as much as it did in the Penn State model. The SSI 
model did not forecast the troughing that the Penn State 
model did in the northeast corner, although an area of 
lower surface pressure did exist there. The 12-hour fore- 
cast uSing the SSI scheme for KIM and PAMELA developed 
strong southerly flows to the east of the storm at middle 
levels which were similar to those produced by the Penn 
State model using V-bogus 1. However, the SSI model also 
developed strong northerly flows to the west of both storms 
at middle and upper levels. These winds acted to reduce 
the northward movement of the storm. The lack of the 
northerly winds, coupled with the trough in the northeast 
corner and the strong upper level anticyclone, allowed the 
southerly winds to move the storm a greater distance to 
the north in the Penn State model. 

V-bogus 3 and Rh-bogus 2 decreased the upper level 
pmereyclonic flow and strong southerly flow in all cases 
using the Penn State model. Although all cases continued 
to experience convection in the northeast corner and along 
the northern boundary, development of a strong trough was 
reduced and limited to forecast times in excess of 24 
hours. Use of V-bogus 3 and Rh-bogus 2 resulted in improved 
tracks for all cases using the Penn State model. Insertion 
of V-bogus 1 in the SSI model reduced the errors in storm 
track positions for OLIVE and PAMELA and increased the 


errors for KIM (Hodur, personal communication). Rh-bogus l 
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improved the forecast positions in the SSI model for PAMELA 
Slightly over those using V-bogus 1 without moisture effects. 

The Penn State model produced smaller errors than the 
SSI model in storm positions at 12 hours for most of the 
five cases. After 12 hours, the Penn State model using 
V-bogus 3 and Rh-bogus 2 produced smaller errors than the 
SSI using FNWC analyzed winds (no bogus), but larger errors 
than the SSI with V-bogus 1. The exception was Typhoon 
KIM for which all versions of the SSI scheme produced 
smaller errors than the Penn State model after 12 hours. 
Position vector errors are summarized in Table 5. 

Finer resolution was probably responsible for the over- 
all better accuracy obtained with the SSI model. The 
improved resolution allowed the storm to interact with 
small-scale circulations. With less resolution in the 
Penn State model, the large vortex storm was subjected 
primarily to the large-scale steering currents, and 
particularly more steering by the westerlies at upper levels, 
and consequently, more recurvature. 

Both cases of PAMELA produced large errors in later 
forecasts due to insufficient recurvature when V-bogus 3 
Meo muscead fOr anitializgation. The lack of recurvature proba- 
bly resulted from the rapid filling and deterioration of 
Seem Structure which existed in all experiments using 
PAMELA. Observational data indicated that the maximum 


Ss) 
wind speed in the storm was much more intense (60 ms °), 
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yet V-bogus 3 was designed to produce maximum winds of 
22 m = =. AeCoacd mney yey -wOecus 5 was altered to inelinde 


more intense wind speeds by changing equations (8) and (9) 


=O 

vg(ms") = 60 & 0<xr< 60 km (14) 
and 

veims >) = 60 - SE780)30 go <x (15) 
for the same limits of r. 7 bogus (referred to as V- 


bogus 4) assumed the same maximum wind and radius of maximum 
wind as V-bogus 1, but produced a significantly smaller 
een, as Can be seen by comparing Fig. 22 and Fig. 23. 
As a result of a greater vorticity maximum and the smoother, 


V-bogus 4 produced greater maximum winds (67 m a) than 


Weeegus |] (57 m a). However, the greater winds did not 
increase the size of the storm. Because the storm was 

Still reasonably small, boundary values closely resembled 
those of V-bogus 3 and the original analyzed values, whereas 
V-bogus ] had created significant changes in boundary values. 
Resultant vorticities produced by the different bogusing 
techniques are plotted in Fig. 24. Note the faster decline 


in vorticity with increasing radius produced by V-~-bogus 3 


(ama the similar V-bogus 4). 
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48-hour forecasts made uSing V-bogus 4 and Rh-bogus 2 
produced dramatic improvements in the track for Typhoon 
PAMELA (See Fig. 25). By 24 hours, V-bogus 4 produced 
300-400% more convective precipitation and heating over 
the storm, but increased convection in other areas by less 


meen LOZ. 


ee 





v= CONCEUS TONS 


The primitive equation model tested appeared to be 
very sensitive to latent heating effects produced by 
Sonvection. Since the convective parameterization was 
initiated by the vertically integrated moisture conver- 
Gence, proper representation of initial wind and moisture 
conditions at all levels was important. A moisture bogus 
Semoltsting Of a uniform field of low relative humidities 
With an area of higher relative humidities centered over 
the storm enhanced development of the storm while suppressing 
convection and development in areas not associated with the 
Storm. The model also demonstrated some sensitivity to 
small changes in the initial values used in the moisture 
Begus . 

Experiments conducted using the Penn State model and 
enotne= Primitive equation model in conjunction with a 
Simple vorticity bogus (V-bogus 1) showed that forecast 
storm tracks could be improved over forecasts made using 
the FNWC analyzed wind fields. Premature recurvature pro- 
duced in the Penn State model forecasts was retarded by 
uSing a smaller, less intense vortex (V-bogus 3). Alterna- 
tively, more recurvature was induced by increasing the 
Maximum wind used in V-bogus 3. This bogus produced more 
intense winds than V-bogus 1 without introducing a large 


cyclonic area. The smaller more intense vortex improved 
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Peeceast tracks Over V-bogus 1 and the less intense version 
O£ V-bogus 3 in the Penn State model forecasts. It was 
fiemerOre possible to produce significant, and somewhat 
controlled, changes in the forecast tracks of individual 
storms by varying the size and intensity of the initial 
Dogus. 

Miruce experiments should include increased resolution 
of the horizontal domain. The SSI model appeared to have 
an advantage in representing steering currents because of 
Peseriner resolution. The initial wind maximum introduced 
Byeene bogus could not be sustained on the 120 km grid, 
gma marked reductions in intensities occurred for all 
boguses used. Given finer resolution, V-bogus 3 using a 
variable maximum wind based on observations of the storm 
should provide a well-structured storm which does not get 
uao Large. 

The development of strong shears and convergence which 
produced convection along the boundaries, indicated that 
Mo@rEiecactions to the boundary conditions should also be con- 
Sidered in future experiments. Changes to the model which 
could resolve these effects are the inclusion of moisture 
diffusion and the addition of large-scale changes close to 
the boundaries. One possibility 1s the use of a scheme 
described by Perkey and Kreitzberg (1976) where a fraction 
of the large-scale time change (based either on analyses or 


from forecasts produced by a large-scale model) is added to 
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the model predicted tendency depending on the proximity 


Beorene boundary. 
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Beonee le Tyohioon KIM 24-n forecast VF 77111012 
600 mb winds (m s~+) and relative 
humidity (coordinates are grid points 
vice latitude and longitude) 
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150 mb winds (ms ~) and relative 
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Pictte ga. Typhoon KiM initial field SiP-1000 mb 
and winds (ms ~) at o = .95 uSing 
V-bogus 1 (isoline spacing is 4 mb) 
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Figure 3d. Typhoon KIM 36-h forecast SLP-1000 mb 
and winds (m s~t) at o = .95 using 
V-bogus 1 and Rh-bogus l 
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Figure 4. Typhoon OLIVE 36-h forecast VT 78042300 
SLP-1000 mb and winds (m s71) at o = .95 
using V-bogus 1 and Rh-bogus l 
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Figure 5. 
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Typhoon OLIVE 36-h forecast VT 78042300 
accumulated convective precipitation 

(em) using V-bogus 1 and Rh-bogus 1. 
Center value of Olive is 20 cm 


(1soline 
spacing is 1 cm) 
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Figure 6. Typhoon OLIVE 36-h forecast VT 78042300 
convective heating (°C/day) using 
V-bogus 1 and Rh-bogus 1 (isoline 
Siacimicn lamco. e/a) 
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Piaace /. Lyphoon OLIVE 36—-h forecast _ VT 78042300 
SLP-1000 mb and winds (ms -) at 
o = .95 using V-bogus 1 and Rh-bogus 2 
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Figure 8. Typhoon OLIVE 36-h forecast VT 78042300 
accumulated convective precipitation (cm) 
using V=-bogus 1 and Rh-bogus 2 
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Figure 9. 





Typhoon OLIVE 36-h forecast VT 78042300 
convective heating (°C/day) using 
V-bogus 1 and Rh-bogus 2 
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iouimen toa. | Typhoon, KEM E2—h ferecase Vi 77211000 
accumulated convective precipitation (cm) 
using V-bogus 1 and Rh-bogus 1 before 
making changes to model constants (see text) 
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Figure l10b. 
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Typhoon KIM 12-h forecast Vi 77 AAA Oey 6 
accumulated convective precipitation (cm) 
using V-bogus 1 and Rh-bogus l after 
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Figure l10c. Typhoon KIM 12-h forecast VT 77111000 
convective heating (°C/day) before 
changes to model constants (see text) 
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